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THERMAL  DIFFUSION  IN  CARBON/CARBON  COMPOSITES 


I.  INTRODUCTION 

Thermal  conductivity  is  the  property  used  to  relate  heat  fluxes  to 
steady  state  temperature  gradients.  The  application  of  standard  thermal 
conductivity  test  methods  to  heterogeneous  composites  presents  no  parti¬ 
cular  difficulties,  beyond  the  usual  ones  associated  with  this  relatively 
difficult  property  determination. 

The  heat-balance  equation  for  transient  conditions  may  be  written  as 


V  .  XVT  +  internal  sources  and  sinks 


■p  at  * 


where  X  is  the  thermal  conductivity,  C^  is  the  specific  heat  at  constant 
pressure,  and  p  the  density.  If  there  are  no  internal  sources  and  sinks, 

V  •  XVT  «  C  p4?.  (2) 

p  at 

For  homogeneous  materials  whose  thermal  conductivity  is  nearly  independent 
of  temperature,  we  may  treat  X  as  a  constant.  Then  V  •  XVT  becomes  XV2T  and 
equation  (2)  can  be  written  as 

W=T  .  Cp  pis,  D) 

or  V2T  “  Cpp  ■  1  dT 

X  adt’  W 

where  a  =  X/C^p  is  the  thermal  diffuslvity. 


For  one-dimensional  heat  flow 


d2T  dT 


For  heterogeneous  materials,  X  is  not  independent  of  position  and  X 
should  not  be  moved  from  behind  the  Dell  operator  [Eq.  (2)  to  Eq.  (3)].  In 
principle,  then,  the  concept  of  diffuslvity  is  inapplicable  to  heterogeneous 
materials.  However,  in  practice,  materials  are  never  homogeneous  as  point 
defects,  dislocations,  grain  boundaries,  voids,  etc.,  are  present  even  in 
so-called  homogeneous  materials.  Yet  diffuslvity  techniques  have  been  suc¬ 
cessfully  applied  to  these  materials  for  more  than  a  century.  Deliberate 


attempts  to  extend  diffusivity  techniques  to  obviously  heterogeneous  materi¬ 
als  have  been  made  for  many  years.  These  efforts  have  been  intensified 
since  1960's  due  to  the  development  of  the  flash  method  for  measuring  dif¬ 
fusivity  because  of  the  relatively  simple  sample  geometry,  rapidity  of 
measurements,  and  ease  of  handling  materials  with  a  large  range  of  diffusiv¬ 
ity  values  over  large  temperature  intervals  with  a  single  apparatus. 

The  flash  method,  in  which  the  front  face  of  a  small  disc-shaped  sample 
is  subjected  to  a  short  laser  burst  and  the  resulting  rear  face  temperature 
rise  is  recorded,  is  used  in  over  80%  of  the  present  thermal  diffusivity 
measurements  throughout  the  world.  A  highly  developed  apparatus  exists  at 
TPRL  and  we  have  been  involved  in  an  extensive  program  to  evaluate  the  tech¬ 
nique  and  broaden  its  uses.  The  apparatus  consists  of  a  Korad  R2  laser,  a 
high  vacuum  system  including  a  bell  jar  window  for  viewing  the  sample,  a 
tantalum  tube  heater  surrounding  a  sample  holding  assembly,  a  spring-loaded 
thermocouple  and  IR  detector,  appropriate  biasing  circuits,  amplifiers,  A-D 
converters,  crystal  clocks,  and  a  minicomputer-based  digital  data  acquisition 
system  capable  of  accurately  taking  data  in  the  40  microsecond  and  longer 
time  domain.  The  computer  controls  the  experiment,  collects  the  data,  calcu¬ 
lates  the  results,  and  compares  the  raw  data  with  the  theoretical  model. 

The  flash  method  has  now  been  applied  very  successfully  to  composites 
consisting  of  a  second  phase  more-or-less  randomly  dispersed  in  a  continuous 
matrix  [1,2]  and  also  to  composite  materials  composed  of  homogeneous  layers 
sandwiched  together  [3,4]. 

A  third  class  of  composite  materials  is  composed  of  fiber  reinforcements 
imbedded  in  a  more-or-less  homogeneous  matrix.  Such  composites  are  becoming 
Increasingly  important  in  new  technological  applications.  One  major  sub-group 
of  this  class  of  coaiposites  are  known  as  carbon/carbon  materials.  Carbon/carbon 
materials  consist  of  arrays  of  graphite  fibers  lined  up  in  one  direction  im¬ 
bedded  in  a  matrix  consisting  of  arrays  of  graphite  fibers  lined  up  in 
more-or-less  perpendicular  directions  with  the  spaces  filled  with  graphite. 

A  variety  of  geometries  have  been  fabricated  including  1-D,  2-D,  3-D  and 
multidimensional  arrays  depending  on  the  orientation  of  the  fibers.  These 
composites  are  being  used  in  an  increasing  number  of  applications  [5],  at 


least  some  of  which  Involve  transient  heating  conditions.  It  Is  therefore 
of  Interest  to  examine  the  response  of  such  heterogeneous  composites  to  trans¬ 
ient  heat  fluxes . 

The  graphite  fibers  have  high  thermal  conductivity/diffusivity  values 
while  the  matrix  materials  have  relatively  low  values  compared  to  those  of 
the  fibers.  Thus  the  fibers  oriented  in  the  direction  of  heat  flow  act  as 
preferred  paths  for  heat  transfer.  When  such  composites  are  subjected  to 
an  instantaneous  heat  pulse  on  one  surface,  the  temperature  wave  is  not 
planar.  The  results  of  thermal  diffusivity  methods  and  the  applicability  of 
the  concept  of  diffusivity  for  such  materials  legitimately  can  be  questioned. 

A  preliminary  study  [6]  on  a  coarse-weave  3-D  carbon/ carbon  composite  has 
shown  that  the  normalized  rear  face  temperature  response  curve  following  a 
laser  pulse  did  not  follow  the  theoretical  pattern  for  homogeneous  materials. 
However,  the  response  curve  for  a  fine-weave  3-D  carbon /carbon  material  fol¬ 
lowed  the  theoretical  curve  very  closely  [7]. 

A  detailed  experimental  study  of  thermal  diffusion  along  the  principal 

TM 

axes  in  an  AVCO  radially-pierced  fabric  carbon/carbon  material  (HEPN-1) 
has  been  performed  recently  [8].  This  work  was  supplemented  by  computer- 
aided  analysis  of  thermal  diffusion  In  HEPN-1  by  J.  Jortner  (included  as  an 
appendix  in  Reference  8.)  HEPN-1  can  be  considered  to  be  a  2-D  carbon/carbon 
material.  From  the  combined  results  of  the  study  on  HEPN-1  and  the  previous 
research  [6,7],  the  important  parameters  governing  the  results  of  a  particular 
diffusivity  experiment  can  be  delineated.  These  parameters  are  listed  in 
Table  1. 

TABLE  1 

Parameters  Affecting  Thermal  Diffusivity  Experiments 

1.  Magnitude  of  diffusivity  of  fibers  and  of  matrix. 

2.  Fiber  fraction  ratio  in  direction  of  heat  flow. 

3.  Thickness  of  sample  in  relation  to  fiber  bundle  spacing. 

4.  Rear  face  temperature  sampling  area  and  location. 


Furthermore,  from  these  experiments,  it  has  been  shown  that: 


(1)  It  is  possible  to  measure  diffusivity  values  which  correspond  to 
thermal  conductivity  values  by  using  a  sufficiently  thick  sample 
and  large  viewing  area.  The  sample  thickness  required  becomes 
less  as  the  temperature  is  increased. 

(2)  It  is  possible  to  measure  in-situ  diffusivity  values  of  fibers  and 
of  matrix  by  using  point  source  temperature  sensors  and  thin  samples 

(3)  It  is  possible  to  determine  the  effects  of  imperfections  such  as 
delaminations  on  the  thermal  diffusivity.  These  effects  are 
manifest  by  matrix  values  lower  than  those  for  regions  where  the 
defects  are  not  present  (point  source  detector  methods)  and  by  an 
increase  in  diffusivity  values  for  artifacts  measured  in  gaseous 
environments  as  opposed  to  vacuum. 

The  present  study  was  undertaken  to  further  examine  the  use  of  the 
flash  diffusivity  technique  for  coarse-weave  carbon/carbons.  This  study 
consisted  of  measuring  the  diffusivity  of  a  variety  of  samples  machined 
from  on  billet  of  3-D  carbon/carbon  including  measurements  along  the  princi¬ 
pal  axes  and  at  45°  to  the  principal  axes,  intercomparing  the  resultss  and 
examining  the  expected  rear  face  temperature  rise  curves  based  on  finite 
element  techniques. 

II.  SAMPLE  DESCRIPTION 

A  piece  of  carbon/carbon  from  billet  F-ll,  from  the  7-inch  Mantech 
program  [9]  was  used  as  the  sample  material.  The  piece  was  from  Ring  No. 

17  near  the  outside  diameter  in  the  forward  part  of  the  billet  [10].  The 
dimensions  were  5. 8  cm  in  the  axial  direction,  1.8  cm  in  the  radial  direction 
and  2.3  cm  in  the  circumferential  direction.  Thermal  diffusivity  samples 
were  machined  in  the  axial,  (z)  radial  (r)  and  circumferential  (6)  directions. 
In  addition,  samples  were  machined  at  45°  between  the  radial  and  axial 
directions  (rz)  and  at  45°  between  the  axial  and  circumferential  directions 
(z6).  The  nomenclature  and  directions  tested  as  shown  in  Figure  1  (also 
included  in  Appendix  I).  The  sample  dimensions,  masses  and  bulk  density 
values  are  listed  in  Table  2.  Measured  unit  cell  dimensions  and  calculated 
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A.  NOMENCLATURE 


Figure  1.  Nomenclature  for  3-D  Cylindrical  Unit  Cell  and  Sketch 

Showing  Five  Directions  in  Which  Diffusivity  was  Measured. 


TABLE  2 


BULK  DENSITY  VALUES 


SAMPLE 

THICK 

WIDTH 

WIDTH 

MASS 

DENSITY 

(NO.) 

(IN.) 

(IN.) 

(IN.) 

(GMS) 

(GMS  CNT3) 

Axial 

0.3008 

0.4728 

0.5000 

2.2308 

1.914 

Radial 

0.2628 

0.4990 

0.5004 

2.0430 

1.900 

Circumferential 

0.2998 

0.4906 

0.5009 

2.3017 

1.907 

Axial- 

Circumferential 

0.3010 

0.4993* 

— 

1.8672 

1.934 

Axial- 

Radial 

0.2990 

0.5000* 

— 

1.8590 

1.932 

* 


Diameter 
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yarn  bundle  fractions  for  the  principal  ax^  are  given  in  Table  3  (Table  1  of 
Appendix  I).  Thus  the  actual  yarn  bundle  fractions  in  the  axial,  radial  and 
circumferential  directions  for  these  samples  are  0.195,  0.067  and  0.277,  res¬ 
pectively. 

Energy  Materials  Testing  Laboratory  measured  the  thermal  diffusivity  of 
axial  and  radial  samples  from  the  forward  and  aft  sections  of  Rings  17  and  18 
from  billet  F-ll.  These  samples  are  listed  in  Table  4. 

III.  RESULTS 

Thermal  diffusivity  results  were  obtained  on  the  samples  listed  in  Table 
2  over  the  temperature  range  from  about  100  to  1000°C.  For  each  data  set, 
diffusivity  values  were  calculated  at  10,  20,  25,  30,  33  Vs,  40,  50,  60, 

66  V 3 ,  70,  75,  80  and  90%  of  the  rear  face  temperature  rise.  The  rear  face 
temperature  rise  was  one  to  two  degrees  Celsius.  For  homogeneous  materials 
the  diffusivity  values  calculated  at  these  various  percentage  rises  would  all 
fall  within  a  ±2%  band.  However,  for  materials  with  preferred  heat  diffusion 
paths,  the  diffusivity  values  calculated  at  various  percent  rises  decrease 
with  increasing  percent  rise.  The  magnitude  of  this  increase  is  determined 
by  the  parameters  listed  in  Table  1.  In  order  to  indicate  the  magnitude  of 
the  change  in  diffusivity  values  with  increasing  percent  rise,  diffusivity 
values  calculated  at  20,  50  and  80  percent  rise  times  are  given  in  Table  5. 

These  results  are  plotted  in  Figure  2.  From  Table  5  and  Figure  2,  the  following 
observations  are  made: 

1.  The  diffusivity  values  for  the  circumferential,  axial- circumferential 
and  axial-radial  samples  are  relatively  independent  of  percent  rise 
at  all  temperatures. 

2.  The  diffusivity  values  for  the  axial  sample,  especially  those  calcu¬ 
lated  from  50  to  80%,  become  relatively  independent  of  percent  rise 
above  200° C. 

3.  The  diffusivity  values  for  the  radial  sample  decrease  rapidly  with 
increasing  percentage  rise,  especially  at  the  lower  temperatures  and 
between  20  and  50%. 
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TABLE  5 


THERMAL  DIFFUSIVITY  RESULTS 


SAMPLE 

Axial 


Radial 


Circumferential 


Axial- 

Circumferential 


DIFFUSIVITY  VALUE  AT 


TEMP. 

20% 

50% 

80% 

(  C) 

126 

0.690 

0.560 

0.545 

277 

0.503 

0.386 

0.366 

475 

0.375 

0.284 

0.269 

640 

0.333 

0.249 

0.234 

804 

0.289 

0.211 

0.199 

960 

0.191 

0.177 

0.187 

121 

0.530 

0.437 

0.368 

308 

0.344 

0.274 

0.229 

468 

0.261 

0.211 

0.183 

647 

0.216 

0.174 

0.154 

831 

0.178 

0.146 

0.134 

997 

0.159 

0.135 

0.123 

1150 

0.146 

0.126 

0.118 

102 

0.869 

0.830 

0.818 

213 

0.680 

0.624 

0.602 

424 

0.438 

.  0.407 

0.395 

641 

0.339 

0.320 

0.313 

822 

0.287 

0.272 

0.269 

1012 

0.263 

0.242 

0.241 

1198 

0.233 

0.225 

0.230 

103 

0.627 

0.627 

0.612 

233 

0.449 

0.450 

0.449 

367 

0.345 

0.340 

0.330 

552 

0.268 

0.261 

0.261 

691 

0.232 

0.229 

0.228 

97 

0.526 

0.502 

0.507 

165 

0.362 

0.350 

0.330 

321 

0.264 

0.258 

0.251 

495 

0.215 

0.202 

0.199 

656 

0.170 

0.168 

0.175 

730 

0.166 

0.160 

0.163 

Axial-Radial 


0.8 


DIFFUSIVITY  RESULTS 


In  fact,  the  rear  face  temperature  rise  curves  for  the  off-axis 
samples  followed  the  theoretical  models  very  well.  This  is  illustrated  in 
Figure  3  which  shows  the  comparison  between  the  theoretical  model  and  the 
actual  experimental  data  for  the  axial-circumferential  sample  at  552°C. 

Thus  unambiguous  diffusivity  values  are  obtained  for  the  off-axis  samples 
and,  as  will  be  shown  later,  this  will  permit  the  development  of  a  technique 
to  obtain  unambiguous  on-axis  diffusivity  values. 

In  order  to  further  explore  the  diffusivity  values  in  the  axial 
direction,  a  second  sample  was  machined.  This  sample,  here  after  identified 
as  Axial  2  was  0.5  inches  thick  compared  to  only  0.3  inches  thick  for  the 
axial  sample.  The  diffusivity  values  measured  on  Axial  2  are  given  in  Table 
6.  Since  these  results  are  similar  to  those  obtained  for  the  other  axial 
sample,  they  are  not  shown  in  Figure  2.  In  addition  to  the  diffusivity 
values  obtained  using  an  IR  detector,  diffusivity  values  were  measured  on 
Axial-2  using  a  point  source  thermocouple.  The  thermocouple  locations  in 
relation  to  fiber  bundle  positions  are  illustrated  in  Figure  4.  The  diffusi¬ 
vity  values  calculated  at  20,  50  and  80%  rise  times  at  these  positions  are 
given  in  Table  7.  Positions  1,  2  and  6  were  partially  on  axial  fiber 
bundles,  while  3,  4  and  5  were  mainly  on  transverse  fibers.  In  particular 
the  diffusivity  values  for  positions  2  and  6  were  higher  than  those  for  the 
other  positions,  especially  at  lower  percent  rise  times.  At  the  longer  rise 
times,  all  the  diffusivity  values  are  much  closer  together.  In  fact  the 
spread  in  diffusivity  values  is  250%  at  20%  rise  and  only  35%  at  80%  rise. 

The  longer  rise  time  values  agree  well  with  the  value  extrapolated  to  23°C 
from  the  higher- temperature  IR  detector  data  (Tables  2  and  6). 

In  one  experiment,  the  thermocouple  was  well-centered  on  an  axial  fiber 
bundle.  The  results  are  given  in  Table  8.  The  experimental  rise  curve  data 
closely  followed  the  theoretical  model  and  the  diffusivity  values  were  all 
in  the  range  from  2.71  to  2.83  cm2  sec  1  from  20  to  80%  rise.  These  values 
correspond  to  the  in-situ  fiber  bundle  at  23°C. 

The  diffusivity  values  reported  by  EMTL  for  their  axial  and  radial  samples 
are  given  in  Table  9.  These  values  were  calculated  at  50%  rise  times.  The 
EMTL  results  for  the  axial  samples  are  compared  to  the  present  results  in 
Figure  5  and  for  the  radial  samples  in  Figure  6.  For  both  types  of  samples, 
the  present  results  parallel  the  EMTL  results  but  are  substantially  larger. 


FIGURE  3.  NORMALIZED  REAR  FACE  TEMPERATURE  RISE 

FOR  AXIAL-CIRCUMFERENTIAL  SAMPLE  AT  552°C 
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TABLE  6 


THERMAL  DIFFUSIVITY  RESULTS 
(Axial  2) 


TEMP. 

(°C) 

20% 

DIFFUSIVITY 

50% 

RESULTS  AT 
80% 

94 

0.637 

0.614 

0.607 

220 

0.519 

0.437 

0.416 

306 

0.415 

0.345 

0.321 

414 

0.326 

0.296 

0.293 

520 

0.314 

0.276 

0.271 

-Source  Temperature  Sensing  Probe 
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TABLE  7 

DIFFUSIVITY  VALUES  AS  A  FUNCTION  OF  POSITION 
(Axial-2  Sample,  23°C) 


POSITION 

CALCULATED 

DIFFUSIVITY  AT 

(NO.) 

20% 

50% 

80% 

1 

1.418 

1.042 

0.925 

2 

2.252 

1.678 

1.153 

3 

1.210 

1.005 

0.958 

4 

0.879 

0.775 

0.782 

5 

1.017 

0.892 

0.828 

6 

2.104 

1.401 

0.978 

1 


TABLE  8 


THERMAL  DIFFUSIVITY  VALUES  FOR  A  FIBER  BUNDLE  EXPERIMENT 

THERMAL  DIFFUSIVITY 


23C 

VALUE  TIME (SEC) 

BASELINE:  0.94000  0.  020680 

HALFtfAX:  3.46995  0.020033 

MAXIMUM:  5.99990  0.073100 


TIME (SEC) 
0.  009761. 
0. 012177 
O. 013333 
0.  014610 
0.  015530 
O.  017351 
0.  020033 
O.  023696 
0.  026136 
0.  027239 
O.  029499 
0.  033231 
0.  042620 


FINITE  PULSE  TIME  CORRECTION 


ALPHA 

2.  68307 

2.  74091 
2.  75494 
2.  74287 
2.  72871 
2.  71593 
2.  74435 
2.  71226 
2. 74445 
2.  79052 
2.  82674 
2.  78000 
2.  82117 


PERCENT 
10.  0 
20.  O 
25.  0 
30.  0 
33.3 
40.  0 
50.  0 
60.  0 
66.  7 
70.  0 
75.  0 
BO.  0 
90.  0 


VALUE 
1.  44599 

1.  95198 

2.  20490 
2.  45797 
2.  62663 

2.  96396 

3.  46995 

3.  97594 

4.  33327 

4.  48193 
4. 73493 
4. 98792 

5.  49391 


2.  809385 
2. 776030 
2.  646231. 


25.  0 
50.  0 
75.  O 


TABLE  9 


DIFFUSIVITY  RESULTS  (EMTL) 


EMTL  SAMPLE 
(NO.) 

F-11-F-5-0T-R-18-D-1 


SAMPLE 

(NO.) 


F-  11-F-  5-OB-R-  18-D-  2 


F-  11-F-  5-OT- A- 1 8-I>- 1 


Ax-1 


F-11-F-5-OB-A-18-D-2 


Ax-2 


TEMP. 

DIFFUSE 

(#C) 

(cm  2  s< 

25 

0.515 

49 

0.475 

99 

0.392 

154 

0.320 

253 

0.244 

355 

0.191 

456 

0.161 

553 

0.137 

703 

0.120 

854 

0.104 

1000 

0.094 

24 

0.534 

50 

0.460 

102 

0.359 

151 

0.308 

255 

0.230 

351 

0.183 

451 

0.153 

555 

0.132 

704 

0.110 

851 

0.097 

995 

0.088 

25 

0.653 

55 

0.593 

104 

0.513 

157 

0.441 

256 

0.327 

353 

0.255 

451 

0.210 

555 

0.178 

698 

0.153 

848 

0.136 

988 

0.127 

23 

0.786 

50 

0.682 

99 

0.531 

152 

0.444 

249 

0.346 

349 

0.273 

451 

0.223 

550 

0.196 

647 

0.180 

749 

0.160 

886 

0.143 

992 

0.136 

C  03S/g 


Comparison  of  EMTL  and  TPRL  Diffusivity  Results  for  Radial  Samples 


— -I- 
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IV.  DISCUSSION 

It  has  been  shown  that  the  diffusivity  values  calculated  at  the  higher 
percentage  rise  times  (longer  experimental  times)  approach  or  equal  the 
diffusivity  values  corresponding  to  conductivity  (steady  state)  experiments. 
Using  this  as  the  basis  "best"  diffusivity  values  for  the  five  sets  of 
experiments  are  listed  at  selected  temperatures  in  Table  10.  These  values 
are  plotted  in  Figure  7.  Using  Ei^&tion  (g)  of  Appendix  I,  namely 
kRn  =  l/2(k^  +  kj),  values  for  the  45°  off-axis  directions  were  calculated 
from  the  on-axis  values  of  Table  10.  These  calculated  values  are  compared 
to  the  experimental  values  for  the  two  curves  measured-name ly  axial-circum¬ 
ferential  and  axial- radial  in  Table  11.  In  general,  the  agreement  was 
within  10%.  When  we  consider  the  fact  that  the  off-axis  samples  were  not 
machined  according  to  the  recommended  geometry  (Figure  8,  Appendix  I) ,  all 
the  yarn  bundles  in  the  back  face  viewing  area  were  not  directly  heated  at 
the  front  face  and  this  should  lead  to  lower  experimental  values  for  the 
off-axis  samples.  From  Table  11,  it  is  noted  that  the  experimental  values 
are 'in  fact  lower  than  the  calculated  values  with  only  one  exception 
(axial-radial  at  100°C) . 

"Steady-State"  diffusivity  values  for  100,  200,  500  and  1000°C  are 
plotted  as  a  function  of  fiber  fraction  in  the  direction  of  heat  flow  in 
Figure  8.  The  solid  lines  represent  the  fiber  fractions  given  in  Table  3, 

(Case  I)  and  the  diffusivity  values  of  Table  10  and  Figure  7.  The  dashed 
lines  in  Figure  8  represent  the  situation  in  which  the  average  fiber  fraction 
value  of  6.7%  for  the  radial  samples  is  replaced  by  3.5%  (Case  II). 

Deshpande,  and  Bogaard,  and  Taylor  [11]  have  shown  that  the  fiber  fraction  in 
the  IR  detector  viewing  area  may  be  substantially  different  than  that  for 
the  over  all  sample  average,  especially  for  radial  samples.  Since  the  actual 
viewing  area  of  our  IR  detector  is  not  known  accurately,  we  wished  to  examine 
the  effects  caused  by  the  uncertainty  of  the  actual  fiber  fraction  in  the 
viewing  area.  The  diffusivity  values  obtained  by  extrapolating  the  solid 
lines  [Case  I]  and  dashed  lines  [Case  II]  to  0  and  100%  are  plotted  as  a 
function  of  temperature  in  Figures  9  and  10,  respectively.  These  values  corres¬ 
pond  to  those  for  the  matrix  (Figure  9)  and  for  the  fiber  bundle  (Figure  10). 
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TABLE  10 

+ 

"STEADY  STATE"  THERMAL  DIFFUSIVITY  RESULTS 


SAMPLE  DESIGNATION 


TEMP. 

<°C) 

AXIAL 

RADIAL 

CIRCUMFERENTIAL 

AXIAL- 

CIRCUMFERENTIAL 

100 

0.570 

0.352 

0.803 

0.607 

200 

0.408 

0.266 

0.621 

0.479 

300 

0.333 

0.216 

0.483 

0.381 

400 

0.284 

0.182 

0.401 

0.314 

500 

0.252 

0.158 

0.350 

0.275 

600 

0.228 

0.143 

0.315 

0.248 

700 

0.210 

0.131 

0.291 

0.226 

800 

0.196 

0.122 

0.270 

0.211 

900 

0.185 

0.115 

0.255 

0.200 

1000 

0.175 

0.109 

0.242 

0.192 

AXIAL- 

RADIAL 


0.478 

0.313 

0.262 

0.222 

0.195 

0.174 

0.157 

0.147 

0.136 

0.126 


(33SA  W3)  AlIAISfUdlQ  “ldWfcGHl 
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CALCULATED  AND  MEASURED 


TABLE  11 


OFF-AXIS 

DIFFUSIVITY  VALUES 

AXIAL- 

■RADIAL 

RADIAL- C IRCUMFERENTI j 

MEASURED  CALCULATED 

CALCULATED 

0.478 

0.461 

0.578 

0.313 

0.337 

0.444 

0.262 

0.275 

0.350 

0.222 

0.233 

0.292 

0.195 

0.205 

0.254 

0.174 

0.186 

0.229 

0.157 

0.171 

0.211 

0.147 

0.159 

0.196 

0.136 

0.150 

0.185 

0.126 

0.142 

0.176 

"Steady-State"  Diffusivity  Values  Versus  Fiber 
Fraction  in  Direction  of  Measurement. 


Thermal  Diffusivity  of  the  Matrix  as  a  Function  of  Temperature. 


Thermal  Diffuaivity  of  the  Fiber  Bundles 
Function  of  Temperature 


It  was  found  possible  to  remove  most  of  the  axial  fiber  bundles  from 
Sample  Axial  2  and  to  replace  them  with  matrix.  Temperature  rise  curves 
measured  in  the  resulting  sample  tended  to  follow  the  theoretical  model  fairly 
closely  and  a  diffusivity  value  of  0.32  cm2  sec  1  was  obtained  at  110°C  in 
air.  This  is  close  to  the  Case  II  value,  but  is  not  in  good  agreement  with  the 
Case  I  results  (Figure  9).  Attempts  were  also  made  to  determine  the  effects  of 
slight  displacement  of  the  radial  sample  at  100°C.  The  results  for  those 
displacements  were  0.398,  0.429,  0.390  at  the  50%  rise  times  measured  in  air. 
This  variation  is  too  small  to  permit  a  large  error  in  the  fiber-fraction 
assigned  to  the  IR  viewing  area.  However,  two  interesting  observations  were 
made  during  these  tests;  namely  that  the  diffusivity  values  were  noticeably 
lower  and  less  dependent  on  percent  rise  time  than  the  values  recorded  initially 
(Table  5  and  Figure  2) ,  and  the  diffusivity  values  were  noticeably  lower  when 
measured  in  vacuum.  The  diffusivity  values  measured  in  air  and  calculated  at 
80  and  90%  rise  times  were  less  than  10%  below  the  50%  rise  time  values  when 
measured  in  air.  The  radial  sample  had  been  sent  to  J.  Jortner  for  microscopic 
examination  and  this  procedure  Included  sample  polishing.  If  we  substitute 
the  "Steady-State"  diffusivity  values  measured  on  the  polished  sample,  we 
would  still  be  on  the  dashed  line  (Case  II  Figure  8)  at  least  for  the  100°C 
data.  Thus  it  is  possible  that  the  diffusivity  values  measured  on  the  matrix- 
only  sample  agrees  with  the  diffusivity  values  obtained  by  extrapolating  to 
0%  fiber-fraction. 

Reasonable  agreement  was  obtained  between  the  diffusivity  value  measured 
for  a  fiber  bundle  (Table  8)  and  the  diffusivity  value  extrapolated  to  100% 
fiber-fraction  for  either  the  Case  I  or  Case  II  situation  (Figure  10) . 

The  flash  technique  has  been  shown  to  be  a  powerful  tool  for  studying 
transient  heat  flow  in  carbon/carbon  materials.  As  a  result  of  the  present 
work,  the  following  conclusions  can  be  drawn: 

1.  Thermal  diffusivity  values  measured  on  45°  off-axis  samples  yield 
unique  values  which  define  unambiguously  the  "steady-state"  diffusi¬ 
vity  values  along  the  principal  axes. 

2.  Diffusivity  values  measured  along  the  principal  axis,  especially 
below  500°C  in  the  axial  and  radial  directions,  depend  upon  sample 
thickness,  detector  viewing  area,  and  percent  rise  time  chosen  to 
make  the  calculations.  Nevertheless,  it  is  possible  to  obtain 
diffusivity  values  equivalent  to  the  steady  state  case. 


3.  It  is  possible  to  determine  the  diffusivity  of  the  matrix  and 
of  the  fiber  bundles  both  by  direct  measurements  and  by  extra¬ 
polating  diffusivity  values  measured  on  samples  containing  vari¬ 
ous  known  fiber  fractions  in  the  direction  of  measurement  to  0 
and  100%,  respectively. 

4.  It  is  possible  to  determine  the  diffusivity  values  for  the  matrix 
and  fiber  bundles  as  a  function  of  temperature  using  the  extrapo¬ 
lation  method  given  in  3. 

5.  There  is  a  surface  layer  in  carbon/carbon  in  which  inter-constitu¬ 
ent  temperature  gradients  are  significant  and  beyond  which  they 
are  negligible.  The  material  may  be  treated  as  homogeneous  for 
the  purposes  of  heat  conduction  through  parts  which  are  much 
thicker  than  the  heterogeneity  dimension.  However,  within  a 
surface  layer  subjected  to  large  heat  fluxes,  the  analysis  must 
consider  the  heterogeneity. 

Additional  diffusivity  work  on  carbon/carbon  materials  Is  warranted.  Speci¬ 
fically  such  work  should  take  into  account  certain  deficiencies  in  the  present 
study.  These  deficiencies  include: 

(1)  Insufficient  sample  material  for  machining  samples  along  all 
principal  axes  and  all  45°  off-axis  directions 

(2)  Viewing  area  of  IR  detector  and  relation  to  fiber-bundle  locations 
should  be  controlled 

(3)  Off-axis  samples  should  be  configured  as  shown  in  Figure  8  of 
Appendix  I. 

(4)  Effects  of  handling  samples,  especially  polishing,  and  the  effects 
of  gaseous  environment  should  be  investigated. 
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Introduction  and  Summary 

For  coarse  weave  composites,  significant  errors  can  occur  In  Interpreting 
dlffuslvlty  tests  because  of  temperature  differences  among  the  variously 
oriented  yarn  bundles  In  the  test  specimen  (eg.  Reference  1).  One 
manifestation  of  these  errors  Is  that  the  dl ffuslvlties  calculated  from 
the  back face  temperature  traces  In  flash-method  tests  are  significantly 
affected  by  the  temperature  rise-fraction  used  to  reduce  the  data. 

Previous  efforts  to  measure  the  thermal  dlffuslvities  of  3D  carbon-carbon 
composites  have  used  samples  oriented  to  measure  the  principal  conductivi¬ 
ties  -  that  Is,  In  the  directions  parallel  to  each  of  the  three  orthogonal 
sets  of  yarn  bundles  In  the  composite.  This  article  describes  some  tests 
on  a  representative  coarse  weave  composite  in  which  data  also  were  obtained 
In  directions  at  appreciable  angles  to  the  yarns.  The  di ffuslvlties  Inferred 
from  these  off-axis  tests  show  less  dependence  on  rise  fraction,  and  appear, 
therefore,  to  be  more  reliable  Indicators  of  bulk  composite  behavior  than  the 
on-axis  data.  This  finding  leads  to  a  recommendation  that  off-axis  testing 
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be  Included  routinely  In  the  characterization  of  30  composites.  However, 
because  of  the  limited  number  of  tests  conducted  In  the  effort  described, 
further  testing  appears  necessary  to  firmly  establish  the  validity  of 
the  recommendation.  In  particular,  the  dependence  of  the  off-axis  data 
on  specimen  dimensions  and  proportions  needs  exploration. 


Experimental  Results 

A  piece  of  carbon-carbon  billet  F-ll,  from  the  7-Inch  Mantech  program  (Refer¬ 
ence  2),  was  supplied  to  TPRL  for  the  purpose  of  exploring  the  thermal  dlf- 
fusivlty  In  a  coarse-weave  material.  A  schematic  of  the  unit  cell  of  a  3D 
cylindrical  composite  Is  shown  In  Figure  1.  The  piece,  about  22  cubic  centi¬ 
meters  In  volume,  was  taken  from  the  outer  region  of  the  billet  where  the 
yarn  bundle  fractions  In  each  of  the  three  principal  directions  are  approxi¬ 
mately  (Table  1): 

Vr  ■  .067 
va  «  .195 
Vg  -  .277 

Samples  were  excised  at  TPRL  to  measure  diffuslvlty  In  each  of  five  directions 
the  three  "on-axis"  directions,  along  the  radial,  axial,  and  clrc  yarn 
bundles;  and  In  two  "off-axis"  directions,  45-degrees  between  the  axial  and 
radial  directions,  and  45-degrees  between  the  axial  and  clrc  directions. 

These  directions  are  Illustrated  In  the  lower  half  of  Figure  1.  The  purpose 
of  the  45-degrees  off-axis  tests  was  to  provide  data  which  would  check  the 
self-consistency  of  the  diffuslvlty  results,  as  the  off-axis  data  Is,  In 
principle  (Reference  3),  predictable  from  the  on-axis  data. 

The  tests  were  by  the  flash-laser  method,  described  In  Reference  1,  and  the 
data  were  reduced  at  TPRL  from  the  backface  rise  times  corresponding  to  rise- 
fractions  of  20,  50,  and  80  percent.  These  results  are  listed  in  Table  2. 

The  data  have  been  plotted  In  Figures  2,  3,  and  4  for  rise- fractions  of  20,  50 
and  80  percent,  respectively. 


Figure  1.  Nomenclature  for  3D  cylindrical  Unit  Cell  and  Sketch 

Showing  Five  Directions  in  which  Diffusivity  was  Measured. 
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Table  1 

Measured  Unit  Cell  Dimensions  (mm)  and  Calculated  Yarn 
Bundle  Volume  Fractions  for  On- Yarn  Dlffuslvlty  Specimens 
from  O.D.  Region  of  Billet  F-ll 


• 

Radi al 
Sample 

Axl  al 

Sample 

Clrc 

Sample 

rz 

.41 

- 

- 

re 

.81 

- 

- 

zr 

- 

.84 

- 

ze 

- 

1.67 

- 

Cr 

- 

- 

1.07 

Cz 

- 

- 

.97 

- 

2.18 

2.29 

1.63 

- 

1.63 

3.05 

3.30 

- 

Vi 

.067 

.195 

.277 

1)  Dimensions  listed  are  averages  of  several  measurements  made 

on  polished  surface  of  the  specimen  with  a  cal ibrated-eyeplece 
microscope  at  20X  magnification. 

2)  .Yarn  bundle  volume  fractions  calculated  using  relations  of  this 
type: 

r& 

M  ■ 


Vr 
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TABLE  2 

THERMAL  DIFFUSIVITY  RESULTS  (cm^/sec) 


DIFFUSIVITY  VALUE 

AT 

RISE 

SAMPLE 

TEMP. 

ec) 

20Z 

50Z 

80Z 

FRACTION 

effect,  t 

Axial* 

124 

0.627 

0.542 

0.508 

23.4 

126 

0.690 

0.560 

0.545 

26.6 

277 

0.503 

0.386 

0.366 

37.4 

475 

0.375 

0.284 

0.269 

39.4 

640 

0.333 

0.249 

0.234 

42.3 

804 

0.289 

0.211 

0.199 

45.2 

960 

0.191 

0.177 

0.187 

2.1 

Radial 

121 

0.530 

0.437 

0.368 

44.0 

308 

0.344 

0.274 

0.229 

50.2 

468 

0.261 

0.211 

0.183 

42.6 

647 

0.216 

0.174 

0.154 

40.3 

831 

0.178 

0.146 

0.134 

32.8 

997 

0.159 

0.135 

0.123 

29.3 

23.7 

1150 

0.146 

0.126 

0.118 

Circumferential 

102 

0.869 

0.830 

0.818 

6.2 

213 

0.680 

0.624 

0.602 

13.0 

424 

0.438 

0.407 

0.395 

10.9 

641 

0.339 

0.320 

0.313 

8.3 

822 

0.287 

0.272 

0.269 

6.7 

1012 

0.263 

0.242 

0.241 

9.1 

1198 

0.233 

0.225 

0.230 

1  •  3 

Axial- 

Circumferential 

103 

0.627 

0.627 

0.612 

2.5 

233 

0.449 

0.450 

0.449 

0.0 

367 

0.345 

0.340 

0.330 

4.5 

552 

0.268 

0.261 

0.261 

2.6 

0 

691 

0.232 

0.229 

0.228 

1.8 

Axial-Radial 

97 

0.526 

0.502 

0.507 

3.7 

165 

0.362 

0.350 

0.330 

9.7 

" 

321 

0.264 

0.258 

0.251 

5.2 

495 

0.215 

0.202 

0.199 

8.0 

656 

0.170 

0.168 

0.175 

-2.9 

1.8 

730 

0.166 

0.160 

0.163 

Rise  fraction  effect  Is  defined  as  — - - —  x  100  percent, 

“80 

where  and  “gO  are  the  d^us*v*ty  values  at  20* and  80- 
percent  rise  fraction  respectively. 
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20%  Rite 


Billet  F-ll ,  00  Region 
O  Circ 

A  Circ-  AX  (~  45*  ) 
•  Axial 

A  AX- Rad  (~45  ) 

□  Radial 


in>: 


50  %  Rite 

\ 

Billet  F- II ,  OD  Region 

O  Circ 

A  Clrc-AX(~45°) 

•  Axial 

»  \ 

A  Ax- Rad  (~45*) 

A  X 

□  Radial 

TEMPERATURE  ,*C 

Figure  3.  Thermal  Dif fusivities  Inferred  from  50%  Rise-Fraction 


7 


TEMPERATURE  f°C 

Figure  A.  Thermal  Dif fusivities  Inferred  from  80Z  Rise-Fraction 
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Analytical  Relationships  Between  On-Axis  and  Off-Axis  Conductivities 
By  definition,  heat  flux  by  conduction  Is  related  to  thermal  gradients  by: 


C-Q]  -  DO  CT'] 


(a) 


The  matrix  equation  can  be  stated  more  fully  as: 


<11 

“jl 

<k1 


K1j 

Kjj 

*kj 


<1k 

Kjk 

<kk 


(b) 


For  an  orthotropic  material,  when  1,  j,  and  k  are  the  principal  axes  of  the 
material  (eg,  r,  z,  and  6,  In  a  cylindrical  3D  composite): 


(c) 


The  principal  conductivities,  kP,  kz,  and  kg,  can  be  derived  from  on-axis 
tests  diffusivlty  tests  using  the  relationship: 


[K]  *  cpM  (d) 

To  derive  the  conductivity  In  an  off-axis  direction,  n,  from  the  principal 
conductivities,  the  rotation  of  the  thermal  conductivity  tensor  of  equation 
(c)  gives  (Reference  3): 

Knn  '  t2m  Kt  ♦  ^Vj  *  l?nk  h  M 

For  the  simple  case  where  n  Is  In  the  plane  normal  to  the  k  axis  and  4k  is 
the  angle  between  the  n  and  the  1  axes: 

"nn  *  <*1  *  "j)  ‘“^k  4  Kj  (f> 


9 


If  4>k  Is  chosen  as  45  degrees: 


*nn  “  -g(Kl  ♦  Kj)  (fl) 

Thus,  we  see  that  the  conductivities  measured  In  the  off-axis  tests  should 
be  bounded  by  the  on-axis  conductivities,  and,  when  the  angle  Is  45  degrees, 
the  off-axis  conductivity  Is  the  simple  arithmetic  average  of  the  conducti¬ 
vities  along  the  two  neighboring  principal  axes. 

Also,  It  turns  out  that  three  Independent  off-axis  tests  could  be  used  to 
derive  the  full  thermal  conductivity  tensor.  For  example,  three  equations, 
of  the  form  of  equation  (g),  can  be  solved  simultaneously  to  yield: 


Kr  - 

Krz 

+  - 

*0Z 

*z  " 

Kp^ 

-  Kre  + 

*02 

(h) 

Ke  * 

-Kp| 

+  Kr©  + 

*02 

where  Is  the  conductivity  In  the  direction  45  degrees  between  the  r  and 
z  yarns,  and  Kr@  and  Kgz  are  similarly  defined. 

Discussion  of  Results 

The  measurements  In  the  directions  of  the  yarns  rank  themselves  In  the  expected 
order.'  That  Is,  the  dlffuslvlty  Increases  with  Increased  yarn  volume  fraction. 
Figure  5,  a  cross-plot  of  the  Information  In  Figure  3,  shows  that  the  rise  In 
dlffuslvlty  with  fiber  volume  fraction  Is  approximately  linear.  A  linear  re¬ 
lationship  between  thermal  conductivity  and  yarn  volume  fraction: 

Ki  -AVi  +  B 


has  been  found  to  give  good  correlation  with  a  wide  range  of  data  for  3D  carbon 
carbon  composites  (Reference  4). 
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Rite  Fraction 


J— 20% 

—  50% 
-—80% 


T  / 2.00*0  / 

i  / 

/  /  . 


/£  '/500»C^ 

T  5X 

/  ^  ^900mC 


YA' 

h'' 


YARN  BUNDLE  VOLUME  FRACTION,  Vj 

Figure  5.  Trends  of  Diffusivity  Vs.  Volume  Fraction 
of  the  Yarn  Bundles  in  the  Direction  of 
Measurement. 


42 


Inspection  of  Table  2  shows  that  the  effect  of  rise  fraction  on  measured 
dlffuslvlty  Is  greatest  In  the  radial  and  axial  directions  and  least  In  the 
45-degrees  (off-axis)  specimens.  The  reason  probably  Is  that  the  texture  of 
the  45-degree  surfaces  Is  more  "uniform"  than  the  textures  of  the  on-axis 
specimens.  That  Is,  more  of  the  45-degree  surface  represents  yarn  bundle 
sections  (Figure  6);  thus,  the  mean  distance  between  high  conductivity  areas 
Is  smallest  on  the  45-degree  surfaces,  Implying  that  the  temperature  gradients 
within  these  surfaces  are  smaller  than  within  the  on-axis  surfaces. 

This  effect  might  be  quantified  In  terms  of  the  area  fraction  of  yarn  bundles 
at  the  front  and  backfaces  of  the  specimen.  Figure  7  shows  how  the  rise- 
fraction  effect  (listed  In  the  last  column  of  Table  2)  varies  with  yarn  bundle 
area  fraction.  The  trend  line  has  been  drawn  with  the  aid  of  the  two  extreme 
points  (at  area  fractions  of  zero  and  unity)  that  represent  "homogeneous 
composites"  for  which  the  rise-fraction  effect  should  be  close  to  zero. 

The  trend  suggests  that  using  three  off-axis  specimens.  Instead  of  three 
on-axis  specimens,  would  provide  more  reliable  estimates  of  bulk  composite 
diffuslvltles. 


As  pointed  out  In  the  preceding  section,  complete  characterization  of  the 
conductivity  tensor  of  an  orthogonally  reinforced  composite  requires  only 
three  Independently-oriented  measurements.  From  these,  the  conductivity  In 
any  direction  may  be  calculated.  Thus  the  45-degree  measurements  are.  In 
principle,  redundant.  For  a  coarsely-textured  composite,  however,  the  45- 
degree  measurements  (assuming  they  are  more  reliable  than  the  on-axis  measure¬ 
ments,  as  Implied  In  the  preceding  paragraph)  can  serve  as  a  check  on  the 
consistency  of  the  data,  via  equation  (g). 

Comparing  the  off-axis  results  to  the  on-axis  results,  we  see  that  data  from 
the  shorter  rise-fractions  (Figures  2  and  3)  obviously  are  not  self-consistent 
because  the  45-degree  diffuslvltles,  are  outside  the  ranges  defined  by  the 
two  relevant  on-axis  diffuslvltles.  For  the  80-percent  rise  fraction,  the 
45-degree  results  are  between  the  relevant  on-axis  results,  but  somewhat 


12 


Comparison  of  Textures  in  the  Off-Axis  45°  9Z  Specimen  and  the 
on-Axis  0-Specimen. 
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lower  than  the  mean  value  expected  If  the  data  were  self-cons 1 stent'  and  If  the 
samples  were  properly  cut  at  45-degrees. 

It  should  be  noted  that  some  data  scatter  will  result  from  slight  variations 
In  unit-cell  dimensions  and  yarn  bundle  volume  fractions  from  specimen  to 
specimen;  for  example.  Table  1  shows  variations  of  about  5%  In  unit  cell  di¬ 
mensions  Sr  and  &&  for  the  specimens  used  here.  Thus,  perfect  consistency 
Is  not  to  be  expected  for  this  variable  material. 

The  relatively  low  values  of  the  45-degree  off-axis  dlffuslvltles  may  be  due. 

In  part,  to  Improper  selection  of  the  specimen  proportions.  Because  the 
yarn  bundles  are  the  primary  carriers  of  heat  In  the  composite.  It  would  be 
better  to  proportion  the  specimen  as  Illustrated  In  Figure  8  so  that  all  the 
yarns  In  the  back  face  viewing  area  would  be  directly  heated  at  the  front 
face.  Inspection  of  Figure  6  shows  that  such  was  not  the  case  In  the  tests 
performed  here.  Because  some  of  the  yarns  In  the  actual  viewing  area  were 
heated  only  by  conduction  within  the  composite.  It  Is  likely  that  the  resulting 
value  of  dlffuslvlty  Is  lower  than  would  be  obtained  in  a  wider  specimen. 

Conclusions  and  Recommendations 


The  results  described  above  support  the  conclusion  that  diffuslvltles  derived 
from  off-axis  tests  are  less  subject  to  errors,  arising  from  temperature  grad¬ 
ients  within  a  coarse-textured  composite  specimen,  than  dlffuslvltles  derived 
from  on-axis  tests. 

As  off-axis  tests  can  be  used  to  characterize  completely  the  conductivities  of 
a  30  composite,  we  recommend  that: 

1)  Off-axis  tests  be  used.  Instead. of  on-axis  tests,  to  characterize 
the  conductivities  of  coarse-textured  3D  composites. 

2)  The  off-axis  specimens  be  proportioned  as  shown  In  Figure  8. 
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CRITERION:  d  >  d  +  2t 
e  v 

d  =  diameter  of  area  exposed 
e  to  laser 

d  ■  diameter  of  area  of 
v  temperature  measurement 


Figure  8.  Recommended  Off-Axis  Specimen  Geometry 


3)  The  actual  angles  at  which  the  specimens  are  cut  should  be  measured 
and  equation  (e)  be  used  to  Infer  the  full  conductivity  tensor. 
Instead  of  the  simpler  equation  (g)  which  was  used  here  for  Illus¬ 
trative  purposes. 

Before  the.  preceding  recommendations  are  adopted.  It  would  be  prudent  to  do 
some  more  exploratory  testing  with  emphasis  on  the  effects  of  specimen  geome 
try  (Figure  8)  on  off-axis  data. 
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Nomenclature 


temperature-dependent  constants 
specific  heat 
orthogonal  coordinates 

cosine  of  the  angle  between  the  1  and  n  vectors 

thermal  conductivity  In  the  1-dl recti on,  a  principal  axis 
of  the  material 

thermal  conductivity  constant  relating  heat  flow  In  the 
1-dlrectlon  to  a  temperature  gradient  In  the  j-dlrectlon 

heat  flux  In  the  1-dlrectlon 

radial,  circumferential,  and  axial  directions  In  a  cyllndrlcally 
reinforced  billet. 

temperature  gradient  In  the  1-dlrectlon, 
volume  fraction  of  yarn  bundles  parallel  to  the  1-dlrectlon 
distance  along  the  1-dlrectlon 
thermal  dlffuslvltles,  analogous  to  Kj  and  Kjj 
angle  of  rotation  of  coordinate  system  about  the  k  axis 
mass  density 


ANALYSIS  OF  TRANSIENT  TEMPERATURE  RESPONSE  OF  A  CARBON-CARBON  COMPOSITE 
DURING  CONTINUOUS  HEATING  AT  ONE  SURFACE* 
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February  1982 

INTRODUCTION 


The  study  reported  earlier  (Reference  1),  on  transient  temperature  responses 
of  carbon-carbon  composites  during  flash-di ffuslvity  tests,  has  been  continued 
to  explore  the  effects  of  the  composite's  heterogeneity  during  situations 
In  which  the  heating  Is  continuous,  as  In  a  rocket  motor  firing. 

In  the  previous  study  it  was  found  that  temperature  gradients  in  thermal  dlf- 
fusivity  specimens  of  carbon-carbon  material  can  affect  the  diffusivity  value 
derived  from  a  flash  test,  and  that  increasing  the  axial  length  of  the  speci¬ 
men  decreases  the  errors  that  might  occur  from  this  cause.  Although  the 
average  behavior  of  relatively  large  pieces  of  carbon-carbon  was  shown  equi¬ 
valent  to  the  behavior  of  a  homogeneous  material  having  the  same  average 
diffusivity,  concern  was  expressed  that  the  Interconstituent  temperature 
gradients  close  to  a  heated  surface  may  be  significant  to  localized  thermal 
processes  such  as  ablation.  This  concern  motivates  the  calculations  reported 
here.  That  is,  we  ask  how  significant  are  the  temperature  gradients  in  a 
carbon-carbon  composite  that  Is  heated  continuously  at  one  surface  as  in  the 
case  of  a  rocket  motor  throat  Insert. 
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CASE  ANALYZED 


The  same  "TASC"  computer  program,  described  and  used  In  the 
(Reference  1),  was  applied  to  the  same  sort  of  ax 1  symmetric 
the  composite  (Figure  1),  using  the  same  properties  for  the 
(Table  1).  The  case  analyzed  Is  schematically  described  In 
mesh  used  In  the  analysis  Is  shown  In  Figure  3. 


previous  study 
Idealization  of 
two  constituents 
Figure  2;  the 


The  model  represents  a  single  yarn  bundle  and  Its  proportionate  share  of  sur¬ 
rounding  material.  The  center-to-center  spacing  of  the  yarns  In  the  composite 
Is  0.24  cm  and  the  diameter  of  the  yarn  bundles  Is  0.12  cm.  The  yarn  bundle 
has  a  high  conductivity  along  Its  length  (the  "axial"  direction)  and  a  low 
conductivity  transverse  to  the  length.  The  surrounding  material  has  an  Inter¬ 
mediate  (rather  low)  conductivity  and  Is  assumed  Isotropic.  The  composite 
Is  heated  at  the  front  face  (at  the  top  of  the  sketch  In  Figure  2)  by  radia¬ 
tion  from  a  3000°K  source;  the  radiation  exchange  Is  assumed  to  occur  with 
a  view  factor  of  unity  and  both  the  source  and  the  composite  have  emmlsslvl- 
tles  of  one.  The  composite  Is  2.5  cm  thick  (approximately  one  Inch)  and  the 
back  face  radiates  to  a  300°K  environment  also  with  an  emlsslvlty  of  one. 
Initially,  at  time  equals  zero,  the  composite  Is  at  a  uniform  temperature  of 
300°K.  Heating  begins  when  the  front  face  Is  suddenly  exposed  to  the  3000°K 
source. 


The  use  of  radiative,  rather  than  convective,  heat  Input  was  merely  a  matter 
of  convenience  In  using  the  TASC  code.  Although  convection  Is  probably 
the  dominant  mode  In  a  rocket  motor,  we' felt  that  the  major  aspects  of  Inter¬ 
constituent  temperature  differences  could  be  Illustrated  with  the  case 
analyzed.  The  numerical  value  of  the  coldwall  heat  flux  (that  Is,  the  heat 
flux  to  the  front  surface  at  the  start  of  heatlnq)  Is  2.27  x  10?  W/m2,  which 
Is  of  the  order  of  the  heating  rate  In  a  rocket  nozzle. 

As  before  (Reference  1),  It  should  be  noted  that  the  properties  and  the 
dimensions  used  to  describe  the  composite  are  Intended  to  be  Illustrative  of 
the  general  nature  of  such  materials  and  not  necessarily  specific  to  any 
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3000°  K  SOURCE 


Figure  2.  Continuous  Heating  Case 
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BACK  FACE 


0.02  cm 


Mesh  Used  in  TASC  Calculations.  Circles  shew  Nodal 
Points  for  which  Date  is  Plotted  in  Figure  4. 
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particular  carbon-carbon  composite.  In  keeping  with  the  illustrative  purpose, 
and  to  avoid  unnecessary  effort  in  this  context,  the  materials  properties 
were  taken  to  be  independent  of  temperature. 


RESULTS 


The  results  are  summarized  in  Table  2  and  Figures  4,  5,  6,  and  7. 

Figure  4  shows  the  temperature  histories  of  four  points  in  the  model, 
corresponding  to  the  centers  of  elements  marked  with  small  circles  in 
Figure  3.  In  other  words,  these  points  represent 

a)  the  front  face  near  the  centerline  of  a  yarn  bundle 

b)  the  front  face  near  the  outermost  region  of  the  surrounding  material 

c)  the  back  face  near  the  yarn  bundle  centerline 

d)  the  back  face  near  the  outermost  region  of  the  surrounding  material 

We  see  that  a  significant  temperature  gradient  exists  between  constituents 
at  the  front  face,  throughout  the  period  analyzed,  whereas  the  gradient 
at  the  backface  appears  negligible  (on  the  scale  of  the  plot). 

The  remaining  illustrations  are  Intended  to  show  how  deep  a  region  near  the 
front  face  is  significantly  affected  by  interconstituent  temperature 
gradients. 

Figure  5  shows  the  calculated  temperature  gradients  along  the  specimen  length, 
near  the  yarn  bundle  centerline  and  near  the  outer  edge  of  the  surrounding 
material ,  at  various  times  during  heating.  It  appears  that  the  depth  to  which 
a  "significant"  gradient  exists  Is  about  0.15  cm  at  all  times  plotted. 

Figures  6  and  7  are  provided  to  show  the  Interconstituent  temperature  differences 
directly. 
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Table  2.  Selected  Results  of  Computer  Calculation 


322.2  466.7  784.7 


()1 


Figures  6  and  7  show  that,  close  to  the  front  surface,  the  less  conductive 
constituent  (the  surrounding  material)  Is  hotter  than  the  more  conductive  con¬ 
stituent  (the  yarn  bundle),  whereas  the  situation  Is  reversed  beyond  a  depth 
of  about  0.3  cm.  Also,  the  magnitude  of  the  Interconstituent  gradient  Is  well 
over  100°C  at  the  surface,  but  drops  to  less  than  40°C  more  than  about  0.2 
cm  Into  the  body.  At  depths  greater  than  about  1  cm,  the  gradient  Is  always 
less  than  about  10°C. 

It  appears  that  Interconstituent  gradients  would  be  of  greatest  concern  within 
a  distance  of  the  heated  surface  of  about  a  unit-cell  dimension  (about  0.24 
cm,  the  diameter  of  the  model  analyzed.  In  the  present  case).  As  temperature 
changes  of  the  order  of  100°C  can  be  expected  to  have  measurable  effects  on 
ablation,  we  conclude  that  differences  In  the  conductivity  of  the  consti¬ 
tuents  are  important  to  consider  In  the  analysis  of  ablation  rates.  In  this 
regard.  It  should  be  pointed  out  that  the  analysis  does  not  directly  address 
ablation,  and  that  it  does  not  consider  the  effects  of  possible  differences 
In  other  relevant  properties  of  the  constituents,  such  as  emlssivlty,  kine¬ 
tics  constants,  internal  surface  area,  etc. 


CONCLUDING  REMARKS 


The  illustrative  analysis  shows  that  there  Is  a  surface  layer  within  vftich 
interconstituent  temperaure  gradients  are  significant  and  beyond  which  they 
are  negligible.  The  particular  case  analyzed  supports  the  intuitive  notion 
(analogous  to  St.  Venant's  principle  In  stress  analysis)  that  the  dimension 
of  the  surface  layer  Is  of  the  order  of  the  dimension  of  the  hetergenelty. 

That  Is,  the  affected  surface  layer  In  this  case  Is  approximately  as  thick 
as  the  center-to-center  distance  between  yarn  bundles.  However,  It  appears 
dangerous  to  generalize  to  this  extent  on  the  basis  of  one  case;  surely  the 
ratio  of  constituent  dimensions,  the  ratio  of  constituent  axial  conductivities, 
the  ratios  of  transverse  to  axial  conductivities,  among  other  factors,  will 
Influence  the  relative  thickness  of  the  surface  layer.  Nevertheless,  the 
results  support  the  view  that  the  material  may  be  treated  as  homogeneous  for 
the  purposes  of  heat  conduction  through  parts  of  thickness  much  greater  than 
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the  heteroqenelty  dimension,  but  that  consideration  of  heterogeneity  Is 
necessary  for  analysis  of  events  within  a  layer  close  to  any  surface  at  which 
heat  Is  transfered  at  high  rates. 
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